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Abstract 



In this paper we derive a number of exact solutions of the discrete equation 
^ ^ .^i^ ,^ . _ -2z^xl + (r? - 3(^-^ - zl)xl + 

Xn+l-^n-l Xn\Xn+l -i- Xn-l) — z ■ ; 77 ; ^ , 

[Xn + Zn+ l){Xn + - J) 

where Zn = n5 and r], 5, fj, and 7 are constants. In an appropriate limit (1) reduces to the fourth 
Painleve (PIV) equation 

(Pw 1 / dwX^ •? . 9 , 9 N /3 

i^ = 2^[al) +l»'+W^+2(.^-„)» + ^, (2) 

where a and P are constants and (1) is commonly referred to as the discretised fourth Painleve 
equation. A suitable factorisation of (1) facilitates the identification of a number of solutions which 
take the form of ratios of two polynomials in the variable Zn. Limits of these solutions yield rational 
solutions of PIV (2). It is also known that there exist exact solutions of PIV (2) that are expressible 
in terms of the complementary error function and in this article we show that a discrete analogue 
of this function can be obtained by analysis of (1). 
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1. Introduction 

The discrete fourth Painleve equation (d-PIV) was recently derived by Ramani, Grammaticos 
& Hietarinta [1] using the method of singularity confinement [2]. It is given by the three-point, 
non-autonomous mapping 

, . , ^ _ -2znxl + (r? - 3S-^ - zl)xl + jj? 

with the variable x„ to be found in terms of Zn = n5 + C,. This mapping is identified as the 
discretised version of the continuous fourth Painleve equation (PIV) 

d^W 1 { ,33, . 2 , o/_2 . , 



+l«'^ + W + 2(.^ -„)„ + -, (2) 

where a and /3 are arbitrary constants. This is observed by taking the limit of (1) as (5 — 0, with 
7 = 1/5 and and /x finite. This process yields (2) with the parameters a and (3 in that equation 
related to rj and ji according to 

Oi = \ri and (3 = (3a,b) 

The six (continuous) Painleve equations (PI-PVI) were first derived around the turn of the 
century in an investigation by Painleve and his colleagues into which second-order ordinary 
differential equations have the property that the singularities other than poles of any of the solutions 
are independent of the particular solution and so are dependent only upon the equation (cf., [3]); this 
property is now known as the Painleve property. There has been considerable interest in Painleve 
equations over the last few years primarily due to the fact that they arise as reductions of soliton 
equations solvable by inverse scattering as first demonstrated by Ablowitz and Segur [4] . Although 
first discovered from strictly mathematical considerations, the Painleve equations have appeared 
in various of physical applications (ef., [5] and the references therein). The Painleve equations 
may also be thought of as nonlinear analogues of the classical special functions though they are 
known to be transcendental since their solution is not expressible in terms of elementary functions. 
However rational solutions and one-parameter families of solutions of the Painleve equations are 
expressible in terms of special functions are known to exist for particular values of the parameters. 
For example, there exist solutions of PIT, PHI and PIV that are expressed in terms of Airy, Bessel 
and parabolic cylinder functions, respectively (cf., [6]). 

Recently there has been considerable interest in integrable mappings and discrete systems, 
including discrete analogues of the Painleve equations. Some of these mapping and discrete 
equations arise in physical applications. For example, a discrete analogue of PI (d-PI) arose in the 
study of the partition function in a two-dimensional model of quantum gravity [7-10]. Subsequently 
a discrete analogue of PII (d-PII) was derived in [11,12] and later discrete analogues of PIII-PV (d- 
Plll-d-PV) were obtained [1]; for further details on the derivation of the discrete Painleve equations 
see, for example, [13,14]. One important result of the investigations is that the form of the discrete 
Painleve equations is not unique since there exist several possible discrete analogues of the Painleve 
equations. Kajiwara et al. [15] and Grammaticos et ah [16] have derived exact solutions of d-PII 
and d-PIII in terms of discrete Airy and discrete Bessel functions, respectively, in analogue to the 
aforementioned results for the associated continuous Painleve equations. We further remark that 
Lax pairs and isomonodromic deformation problems are known to exist for d-PI [10], d-PII [12] 
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and d-PIII [17]. However, at present, there is no discrete analogue of PVI, nor are there Lax pairs 
for the versions of d-PIV and d-PV derived in [16]. 

In recent work [18,19], wc have been concerned with the investigation of Backlund 
transformations and exact solutions for PIV (2) together with an examination of various applications 
of these solutions to several physically motivated nonlinear partial differential equations. The 
purpose of the present article is to discuss some new solutions of d-PIV (1). In [18] we demonstrated 
how all known exact solutions of (2) can be categorised into one of three families; in two of these 
solutions can be determined in terms of the complementary error and parabolic cylinder functions 
whilst the third family consists of solutions which can be expressed as the ratio of two polynomials 
in z. Since d-PIV reduces to (2) in the appropriate limit then it can be expected that exact 
solutions of (1) exist which should tend to known continuous solutions in the same limit. In §2 we 
shall principally concern ourselves with rational solutions of d-PIV and in §3 we identify solutions 
of (1) which can be thought of as discrete analogues of the complementary error function hierarchy 
discussed in [18]. In §4 we make a few closing remarks. 



2. Rational solutions of d-PIV 

Before we start our study of rational solutions of (1) it is convenient to summarise the state 
of knowledge of rational solutions for PIV (2) . It is easily verified that w = Xj z satisfies (2) with 
the parameter choices a = 2, /3 = — 2 and this is a simple example of a rational solution of this 
equation. Two families of such solutions take the forms 

^z;a,/?) = ^^^, ^^;a,/3)=_2^+^^zi(£), (4a,b) 

where Pm{z) and Qm{z) denote some polynomials of degree m consisting of either entirely even or 
else entirely odd powers of z. Murata [20] has shown that (2) admits unique rational solutions of 
type (4a) or (4b) if the parameters are of the form 

{a,l3) = {±k,-2{l + 2n + kf); k,neZ, n < -1, k > -2n, (5a) 

and 

{a,P) = {k,-2{l + 2n + kf); k,neZ, n>0, k > -n, (5b) 
respectively. A further family of rational solutions of PIV is characterised by 



w{z;a,P) = -lz+^^, (6) 



where {a, (3) = (?t-i, — 1(1 + 3n2)^) with ni and n2 either both even or both odd integers. An 
extended discussion of the important properties of these rational solution families together with 
tables containing the first few solutions in each of these three hierarchies are to be found in [18]. 
Tamizhmani, Grammaticos & Ramani [21] noted that if we set 

ri = 35-^ + 7^ - 2{a^ + b^), = a"^ - b^, 

then the d-PIV equation (1) can be factorised as 

{xn + a + b){xn + a-b){xn- a + b){xn- a-b) 

{Xn+l + Xn){Xn-l + Xn) = —^7 7 , (7) 

[Xn + Zn+ jjiXn + Zn - l) 
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and the parameters a. and /?, as defined by (3), are now given by 

« = i^-' + + ^\ ii = -W~ b'f. (8) 

Simple rational solutions of (7) can be found with Xn proportional to Zn', these are 

xn = -2zn, a = i,5 + 7, 6=i(5-7, a = ^ {6-^ - j^) - (3 = -2j^d\ (9) 

and 

x„ = -|z„, a =15 + 7, 6=|,5-7, a = 1(5-2-72) /3 = -|7252. (10) 

In the limit as 5 — > 0, with 7 = 1/5, these discrete solutions tend to 

w{z;0,-2) = -2z, u;(z;0,-|) = -fz, 

respectively, which are the first members of the PIV hierarchies typified by (4b) and (6). If solutions 
of (7) proportional to 1/zn are sought then it is found that 

^^^_5(5±7)^ a=|5±7, b=^5±-f, a = 5-2-72^275- /3 = -252(5±7)2, (11) 

and so, in the limit as 5 — 0, with 75 = 1, we have x„ — w{z; ±2, —2) = ±l/z and these continuous 
solutions are members of the family described by (4a). 

More complicated rational solutions of d-PIV (1) can be deduced by rewriting (7) as the pair 

^-^^-'=^^''^^^Tt^^' (12a) 

Tamizhmani et al. [21] speculated that such a formalism ought to lead to discrete rational solutions 
though they did not present any such solutions. It is a routine calculation to verify that equations 
(12a) and (12b) are compatible if and only if 

a=i5±7. (13) 

[We remark at this stage that (12) is not the only possibility for the splitting of (7). Easy 
generalisations of (12) include the multiplication of the right hand sides of (12a) and (12b) by 
constants C and 1/C respectively or the taking of the factors in different pairings. However, it can 
be shown that in either of these cases we obtain an incompatible couple of equations so that the 
choice of factorisation (12) is not as specialised and restrictive as it might first appear.] 

If the condition (13) is satisfied then we can seek solutions of (1) by finding solutions of the 
simpler form (12b) (or (12a) would do equally well). If we write Xn = Pn/Qn then (12b) can be 
recast as 

Kn+iP„+i = (z„ ± 7 + 5)P„ - fiQn, (14a) 

Kn+lQn+l = -{Zn =F l)Qn " Pn, (14b) 
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where the 'separation' parameter Kn+i can depend both on n and Zn (recall that = — b"^). 
However, if for the moment we take Kn+i = 1) then eliminating P„ between (14a) and (14b), setting 
J = 1/6 and taking the limit as 6 shows that Qn q{z) where q{z) satisfies where 

= (z2 + /,-l)g. (15a) 

If in this equation we set q{z) = r]{^) with ^ = \/2z then we obtain 



|? = (H' + ^M-i)^, (15b) 



which is the parabolic cylinder equation (cf., [22]). It is well known that when ji = —2n this 
equation admits the solution 

r?(0=He„(Oexp(-if), (16a) 
where He„(^) is the Hermite polynomial of degree n given by 

He„(0 = (-l)"exp (if) — [exp (-if)] . (16b) 

If the variable Qn is eliminated from equations (14) and the usual limit taken then it follows 
in a manner similar to that outlined above that rational solutions Xn = Pn/Qn of (1) exist which 
tend to a function of the form — He^(^)/Hem(^) whenever ji = —2m, with m G N. It has been long 
established that solutions of (continuous) PIV (2) can be related to Hermite functions whenever 
the parameters a and P take certain values. In particular, Lukashevich [23] proved that 

w{z;-{u + l),-2i.'') = -^^, (17a) 

(pu az 

w{z- -V, -2{v + \f) = -2z - i-^, (17b) 

w{z;v,-2{y-r\f) = -2z-^^, (17c) 

.(.;. + l,-2.^)^i-^, (17d) 
where (\>v(z) and are any solutions of the Weber-Hermite equations 

2z-p + 2z/^^ = 0, (18a) 



dz^ 

—-^ + 2z— 2vil)^ = 0. (18b) 

dz^ dz 

li v = n, with n a positive integer, then 0„(z) and ipn{z) are polynomials of degree n given by 

d" d" 
^n{z) = (-1)" exp [z^] — [exp {-z^)] , ijn{z) = exp {-z'^) — [exp [z^]] . 

(These are also expressible in terms of the Hermite polynomial He„(^) defined above.) 

Motivated by these comments concerning the continuous PIV case, we return to equations (14) 
for the discrete situation. If we write 

(P„, Q„) = {An, Bn) X {-5rv ( ^n-m5^l \ ^ ^^g^ 
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where T{z) denotes the usual Gamma function, and choose /Cn+i = 1; then we obtain the pair 

(z„ - m(5 =F 7)^n+i + (zn + (5 ± 7)^^ = fiB^, (20a) 

{Zn-m6^-/)Bn+l- {Zn^^)Bn= An. (20b) 

For m € N we can find exact solutions of these equations with and Bn taking the forms of 
polynomials in z„, consisting of either only even or only odd powers, and of degrees m — 1 and m 
respectively. The first few solutions in this hierarchy are 

m = 1, Xn = —-^ — —, (21a) 

36(26 ±' ) [24- 6(36 ± J)] 

"^ = '' = zA2zl-Si8S±3,)] ' ^'^^^ 

46(55 ±2j) [2zl-SinS±3,)] 
44 -45(10(5 ±37)4 +352(3(5 ± 7) (45 ±7)' 
55(3(5 ± 7) [44 - 4<5(13(5 ± 37)^ + 3.5^(3(5 ± 7) (4(5 ± 7)] . . 
' z„ [44 -205(45 ±7) + 52(25652 ±12575 + 1572)] ' ^ 

In each case the corresponding value of fi in (20a) is /i = —6m [(m + 1)5 ± 27] and this leads to 
the respective values of parameters a and /?, as defined by (8), given by 

a=|(5-2-72):p(m + l)75-|52[2m(m + l) + l], P = -^m^S^ [{m + 1)6 ±2jf . (22) 

We emphasise at this stage that the discrete solutions (21) are exact and are valid for any 6 and 7. 
We can recover continuous solutions by letting 7 = 1/5 and taking the limit as 5 — which yields 
a = T("i + 1) and P = —2m?. Then from (21) we obtain 

w(z- ^2,-2) =tK (23a) 
^2;T3,-8)=t^^^, (23b) 

-^^■^^'^-''^-^^S^3^ (^^^) 

.(.;:,6,-50) = +4(l^±l^^. (23e) 
^ ' ' ' z(4z4 + 20^2 + 15) ^ ' 

It is then clear that solutions (21) can be thought of as discrete analogues of the (continuous) 
solutions of PIV taking forms given by (17a) and (17d), with v = m. 
If in (14) we let 

(P„,Q„) = (A„,S„) X ^nj. (^£!i_i^!^j ^ (24) 

then we obtain 

(zn -m6± 7)(^„+i - An) = -l-J-Bn, (25a) 
(zr, -m6± 7)-B„+l ± (Zn + ^)Br, = -An. (25b) 
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As previously, exact polynomial solutions of these can be derived for integer m except that now 
An and are of degrees m + 1 and m respectively. Then the first few solutions of this type are 

m = 0, Xn = —2zn, (26a) 

m = l, ,^ = -^A^li^, (26b) 

2zr,[2zl-5{55T3l)] . 

"^ = '' = 2zl-6{26Tl) ' ^'''^ 

, _ 44 - A6{75 T 37)4 + 3<^^(3>5 ^ ^){26 ^ 7) 

[24 - <^(85 T 37)] ' ^'^^^ 

2z„ [44 - 20(5(3,5 T 7)4 + ^^"(146^^ T 957(5 + 157^)] 
' 44 - 4(5(10(5 T 37)4 + 3,52 (4(5 T 7) (3-5 ^ 7) ' 

The corresponding values of the parameters a and /? are given by 

a= |(,5-2- 7^) ±7717,5 - 4542m(m + l) + l], /? = -i(m + 1)^,5^(777,5 T 27)^ (27) 

and in the limit 6 —> with 7 = l/,5, solutions (26) reduce to the functions 

w;(z;0,-2) = -2z, (28a) 
w{z;±l,-8) = -2z=f^, (28b) 

w{z; ±2, -18) = -2z T (28c) 

.(.;±3,-32) = -2.^|g|l|, (28d) 

M.;±4, -50)^-2.^ J;f;f^3 . (28e) 

These are all members of the so-called '—22;' hierarchy of rational solutions for PIV (2) and are of 
the form (17b) or (17c), with u = m. It is straightforward to obtain further exact rational solutions 
of d-PIV by solving the pairs (20) or (25) for higher values of m. 

We note here that although we have found some rational solutions of d-PIV (1) it can be 
anticipated that further such solutions exist which are not derivable using the procedure described 
above. This deduction follows from the observation that rational solutions of PIV (2) are possible 
for all parameter values as described by (5) and that so far we have discrete solutions which, in 
the appropriate limit, tend to the continuous solutions of type (17). Therefore, there should be 
discrete counterparts to the remaining continuous rational solutions. These can be constructed by 
appealing to some Backlund transformations for d-PIV which were given by Tamizhmani et al. 
[21]. These authors presented a sequence of transformations which, given a solution of (7) with 
parameters a, h and 7, yields a further solution of the same equation but now with parameters 
a + 5,h and 7. In short, suppose that x„ is a solution of d-PIV with parameters a, b and 7. Then 

~ _ XnXn+l + XnjZn + a) + XnjZn - a) + b'^ - , . 

Vn — . ) i^yj 

Xn I Xji-\-l 

with Zn = Zn + also satisfies d-PIV but now with parameters a = a + ^6, 6 = 7 and 7 = 6. The 
subtlety with transformation (29) is that y„ is defined on a lattice of points which is offset by 
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from the original. In order to obtain a solution valid at points coinciding with the initial lattice it is 
therefore necessary to reapply (29) which has the overall effect of raising the value of a by (5 whilst 
keeping b and 7 invariant. Clearly M applications of this sequence will increment a to a + M5 and 
leaves the other two parameters unchanged so that from a starting solution with parameters 

a=|(5±7, 6= (Ar+ i)(5±7, AT = 0,1,2,... (30) 

(which are the parameter values in (1) appropriate to the family of solutions whose first few members 
are listed in (21)) we can obtain a solution with 

a={M+l)S±j, 6= (iV + i)(5±7, (31) 

or, in terms of the parameters a and P as given by (8), 

« = - 7^) ^ 7<^(^ + iV + 1) + [|(M - N){M + iV + 1) - (M + ^f] (5^ (32a) 

P = -|(M - NfS^ [{M + N + 1)S ± 27]^ . (32b) 

In the usual limit our discrete rational solutions will then tend to continuous ones with associated 
parameters a = t(M -j- N + 1), P = — 2(M — N)"^. Now M and N can be chosen so as to force 
these parameters to coincide with the form (5a) for any k and n in the permissible ranges; in this 
way we have a mechanism for deducing discrete analogues of all those rational solutions of PIV (2) 
which take the form Pn-i{z)/Qn{z). 

A similar argument can be applied to the discrete solutions (28) and this demonstrates that 
given these results then suitable application of the Backlund transformations contained in [21] will 
generate another set of exact solutions of (1). These comprise the discrete counterpart to the '—22;' 
rational hierarchy of (2) which itself is characterised by the parameter values given in (5b). 

We remark here that of the three hierarchies of rational solutions for PIV (2), the simple 
procedures outlined above have yielded discrete analogues of only two of these; no solutions 
corresponding to the — |z family (6) have been found. The reason for this is that although the 
full discrete equation (7) admits the solution Xn = —^Zn with a = |5 ± 7, the splitting of this 
equation into the two Ricatti-like forms (12) gives a pair of equations whose compatibility requires 
that a= |<5±7. Thus, as noted by Tamizhmani et al. [21], this solution is not linearizable through 
the splitting assumption and thus it is unsurprising that the discrete analogue of the — |z hierarchy 
of solutions cannot be generated in this way. However, we shall now show how use of the Backlund 
transformation (29) discussed above can be used to increment the values of a or 6 by ±6 and thus 
lead to some more exact solutions in the discrete —^z hierarchy. (It is noted that d-PIV (7) is 
invariant by interchange of parameters a and 6 or by a sign change of either of these. Therefore, 
given that two applications of (29) increases a by (5 it is easy to see how a suitable changes in sign 
of a and b or the swopping of these parameters allows transformations to be found that increase 
or decrease the values of a or 6 by integral multiples of 5.) A few examples of simpler solutions in 
this hierarchy are: 

Xn ~ 3'2n) (33a) 

Xfi = - , (33b) 

oZn 
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^n-- ':f:iz'Txf2 ^ (33c) 



3(2z2 ± 37(5 - 2(52) 

±7-^5, ^ = ^7+1(5, a = f((5-2-72)±27,5-i(52, I3 = -I5\^^5f 



^n-- ...'\^':'-J'ff^f} ... (33d) 



3(24 - 375 - (52) (24 + 375 - (52) ' 



84 - 4^(^ ± 37)4 - 2(52(38(52 t 697(5 + 2772)^ + 9^3(4^3 ^ ^^^^2 ^ ;^q^2^ ^ 3^3) 
"""" 3z„ [44 - 4(5(4(5 ± 37)^2 ± 97(52((5 ± 7)] ' ^ ^'^ 

a = ±7- 15, 6 = ^7 + ^(5, a = |(,5-2-72)±37(5-i(52, /3 = §,52(3(5 27)'. 

In the limit as (5 — > 0, with 7(5 = 1, the solutions (33) reduce to 

wiz-0,-l) = -lz, (34a) 

^z;±l,-|) = -|z±^, (34b) 

^^;±2,-|) = -|^±2i^' (34c) 

= (2.2-3X2.2+3) ' 
..(.•±3 -i) - - 2. ± 3(4z^±4.2 + 3) 

w;(z,±3, gj- 3^ ± ^(4^4 ± 12^2 _ 9)' (34e) 

which belong to the — §2 hierarchy of rational solutions of PIV (2); see Table 4.1.3 of [18] for 
a more extensive list of such solutions. The method of generating the —^Zn discrete rational 
solutions given in (33) is intrinsically less satisfying than that employed for the evaluation of 
solutions (21) and (26) in the other two families; this feature is a direct consequence of the fact 
that the — §.n forms do not arise from a suitable factorisation of d-PIV into a pair of simple 
compatible equations akin to (12). Instead the present method is based on a direct implementation 
of the Backlund transformation (29), which, for the for more complicated solutions in the hierarchy, 
becomes an increasingly laborious task. However, we believe that the Backlund transformation (29) 
will inevitably be needed to make a complete evaluation of each of the hierarchies and provide a 
systematic and efficient procedure for doing this. 



3. Discrete Complementary Error Function Solutions 

In §2 we have deduced several new rational solutions of d-PIV (1). It has been known for a 
considerable time that PIV (2) admits a number of other exact solutions which can be categorised 
into separate families. In addition to the rational hierarchy, solutions can be found which are 
expressible in terms of the complementary error function and, for certain half-integer and integer 
values of the parameters a and in terms of parabolic cylinder functions. We note the close 
connection which exists between the complementary error function and rational solution hierarchies; 
solutions in the former can be found in terms of ^{z) where 

^^/^z) = ^J^, Mz) = 1-C erfc(z) = 1-^ H &M-t^) dt, (35a,b) 

viz) V^i" Jz 
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with C an arbitrary constant. When we choose C = then ^' = and the complementary error 

sohitions reduce to the rational forms whose discrete analogues have been described above. The 
formulation given in §2 enables us to deduce a discretised version of ^{z) and in order to do this 
we return to system (14). If we take = there and select k^+i = —{zn + t) then 

{Zn + l)Pn+l + {Zn-l + S)Pn = 0, (36a) 
{Zn + 7) [Qn+1 - Qn] - Pn = 0. (36b) 



If we choose Zn = nS then 

^n+l — ; 7 

[ 7 + no 
The first of these is easily solved to yield 



Pn and Qn+i - Qn = — (37a,b) 

7 + no 



p - [r(7/<5)]^ p 

- r (-n + 7/<5) r (n + ^/S) ^''^> 

and it is clear that this is an even function in n. The continuous limit of (38) is best obtained by 
setting 'y = 1/S, S = z/n and noting that 

n-l 

(39) 



ln(Pn/Po) = E In (l - ^) - E 1- (l + ^ ) 



Hence it follows that 

lim {ln(P„/Po)} = + O (1/n) , 

n— >oo 

and so 

P„^Poexp(-z2) = p^(2), 

say. Consequently (38) can be thought of as a discrete version of V''(^) ^ exp(— z^) and so the 
solution of (36b) for Q should yield an appropriate discretised form for '(/'(z) which by (35b) is closely 
related to the error function crfc(z). An explicit solution of equation (36b) is not immediately 
apparent, though it can readily be seen that 

n 

^7 + J<^ 

where Qq is a constant and the P^ are as in (38). This sum is not straightforward to evaluate but 
for the present purposes it is adequate to merely note that as 5 — 0, with 7 = 1/5, then 

n 

which in turn asymptotes to the integral of P*(2;); i.e. to a multiple of the complementary error 
function plus a constant. Thus in forming the solution x„ = Pn/Qn we have a discrete equivalent 
to the continuous PIV solution w{z;l,Q) = '^{z) defined above. Solutions in the (continuous) 
complementary error function hierarchy are of the form 

I m ^('■'^(^)) 
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where P and Q are polynomials in z and ^(2;), and exist for precisely those parameter values 
given in (5). What we have developed here is the discrete analogue of the simplest solution in 
the complementary error function hierarchy; furthermore, this is the discrete analogue of the most 
elementary of the "bound-state" solutions of PIV derived in [24] (see also [25]), which themselves 
form a special case complementary error function hierarchy. 

This technique can be adapted to find discrete counterparts to further exact solutions within 
the complementary error function hierarchy. If we take 



4 



Q-v- = ^ (5-2 _ + (m + 1)7(5 - \ [2m(m + 1) + 1] (5^ (41a) 



for m = 0,1,..., and also write the corresponding discretised solution Xn'"-' = P^^/Q^\ then 
equations (14), with the lower sign taken, may be written as 

^Ti^nS = (^n - 7 + 5)Plr^ (42a) 

= -{^n + i)Q^r^ - pt\ (42b) 

where /x'^'") = mS [27 — (m + 1)6]. In the limit as (5 — 0, with = 1, it is clear from (41) that we 
will obtain solutions with a = m + 1, P = —2m? and it is the discrete complementary error function 
solution with index m = that we have derived in (36)-(40) above. If we choose the separation 
variable — ~ i^n + 7 — inS) then by cross-elimination we can show from (42) that the variables 
Pn"*^ and Qn"*^ satisfy 

[z„ + (7 - <5) - m(5]PS - 2(7 - S)Pi";:l - [z^ - (7 - <5) + m(5]P^) = 0, (43a) 
[zn + l-{m + l)(5]Qi^2 - 27Qi7\ - [^n - 7 + (m + mQ^^^ = 0. (43b) 

Now it can be seen that if for some particular choice m = M we write 7 = 7 — (5 in (43a) then 
the governing equation for P^^'' is identical to that for ^-^ with 7 in (43b) written as 7. The 

result of this observation is that given the solution , we can immediately deduce the form 

of Pn^^ simply by replacing all occurrences of 7 in the former solution with 7 — (5. Furthermore, 
since the continuous limit of our discretised solutions is obtained by letting S and j6 ^ 1, 
it is clear that the appropriate limit for Pn^-* is precisely the same as that for Thus, in 

practice, if we have solutions Pi^ and Qn^ then Pji^^ is directly derived from and 
can be obtained by solving (42b). 
The discretised complementary function solutions described here are exact in the sense that 
they are precise solutions of the d-PIV equation (7) but the practical difficulty with them is that 
the derivation of explicit formulae for the solutions is distinctly awkward. The main reason for 
this stems from the fact that a closed form for the Qn^ solution given by (40) is unknown and 
then, because of the nature of systematic manner for deriving further solutions that we have just 
described, this difficulty with the Qn^ form propagates itself through the rest of the hierarchy. 
However, we shall demonstrate how one can obtain the next complementary error function solution 
and this shows how the procedure operates in practice. Using the notation suggested by (41) we 
have the results that 

p(0) _ F i7/^)f p(o) ^(0) _ ^(0) .44, 

-r(-n + 7/(5)r(n + 7/<5) ° ' ^"+1-^° + ^ ^'^^ 
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which we have noted tend to multiples of the functions ip'iz) and il!{z) respectively (see (35)) in 

the appropriate limit. We obtain immediately from (44b) by replacing 7 by 7 — (5, since 

Pn^ = Qn ^ , and so 

n-1 p(0) 

Then we obtain Qn^ by solving (42b) which yields 

»-2 p(l) p(l) 

J 7 



We can readily see the difficulty in simplifying these solutions further but we can show from (45b) 
that as 6 ^ with 7(5 = 1 so 

n-1 

Q^rl^^Sj^P^'^+Q^o^- (46) 

Since Pi^"* ^ ^p{z) in this limit, so (46) indicates that Qn'* asymptotes the integral of ip{z). 
Moreover, since ip{z) satisfies 

d^V r, dV' 



+ 2z-^ = 0, 



dz^ dz 



this then means that 



and hence 



QW^.^(.) + i^ (47) 

(1) 



,(1) = ^ ^ ^ 2 

" -qW 2z^iz)+^'{z) ^{z) + 2z' ^ ' 

where ^(2;) is as defined in (35a). Thus we conclude that Xn^ and Xn^ provide discretised forms 
of the first two complementary error function solutions of the continuous PIV equation (2) with 
parameters as given by (41), i.e. the continuous solutions ^{z) and 2/[^(z) + 2z] (see Table 3.2.2 
in [18]). Further discretised solutions could in principle by generated by proceeding up the chain 
of equations characterised by (42) and (43) although we shall not do this here. Nevertheless, what 
we have shown is how discrete analogues of complementary error function solutions of (2) may be 
obtained. 

The solutions arising from equations (43) are the complementary error function-like 
generalisations of the discrete rational solutions whose first few members are listed in (21). These 
rational forms vanish as z ^ ±00 as do the solutions which are obtained by following the method 
characterised by (41)- (48). However, we demonstrated in §2 how solutions within the so-called 
'— 2z' rational hierarchy can be obtained by appropriate choice of the separation constant in (14); 
now we shall again take /i^'"-' = — {m + 1)5{2^ + mS), m = 0,1, . . . (with the corresponding values 
of a and (3 then as defined in (27)) and the separation constant = -^n + 7 — n^-^- Then, with 
the upper choice of signs, equations (14) show that in this case Pi"^'' and satisfy 

[zn -il + S)- mS] Pj^f^ + 2(7 + <5)Pi7) - [z„ + (7 + <5) + mS] P^) = 0, (49a) 
[z„ - 7 - (m + 1)5]QS2 + 27Q&1 - [^n + 7 + (m + l)<5]Q(r^ = 0, (49b) 

which are identical to (43) under the transformation 7 —7. This then provides the starting 
point for the derivation of further complementary error function-like solutions which are the 
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natural extensions of the '— 22;„' rational forms whose first few members are as in (26). The 
most straightforward example is provided by m = whence, from (44), it can shown that 

n — 1 ^ 

and we have already shown how this function asymptotes to a multiple of tp{z) in the usual limit. 
It then follows that 

{zr. + 7)Pit - (^n + 7 - S)Pi'^ = 27'5Qi°) (51) 
and solution of this equation in the limit as (5 ^ 0, with 7(5 = 1, leads to the result that 

n 

i=o 



Thus Pi^^ asymptotes to the function 



) = -2 y" V(s) ds = -2ztp{z) - 



dV' 
dJ' 



and so the solution 

,(0) 



(o)_Pr , 2zi;{z)+^'{z) 



" g(0) ^^^^ 

and thus we have the first in the sequence of '— 2z' complementary error function solutions of PIV 
(2). A list of the simpler solutions in this hierarchy may be found in Table 3.2.3 of [18] and discrete 
equivalents of these solutions may be derived using the ideas here; the essential ingredient of this 
process is the requirement to solve equations (49) for higher values of m. 



4. Discussion 

The main result of this paper has been the evaluation of a number of exact solutions of the 
d-PIV (1) and consequently answer a number of the questions raised by Tamizhmani et al. [21]. 
These solutions are analogues of known solutions of the continuous PIV (2). This provides further 
evidence of the close relationship and analogy between properties of the continuous and discrete 

Painleve equations. 

We believe that further exact solutions of d-PIV (1) will arise by studying the discrete Backlund 
transformations given in [21]. A detailed study is also required in order to find discrete analogues 
of these solutions of PIV (2) that are expressible in terms of the parabolic cylinder function (z) , 
for non-integer u. An important open question for d-PIV (1) remains that of the derivation of an 
associated Lax; pair and the form of the associated isomonodromic deformation problem, the 
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